Non-locality from N > 2 Independent Single Photon Emitters 



O 
(N 






Oh: 



(N 

00 



in 

o 
o 



^ 



C. Thiel,^ R. Wiegner,i'El J. von Zanthier/ and G. S. Agarwal^ 

^Institut fiir Optik, Information und Photonik, Universitdt Erlangen-Nilrnherg, 91058 Erlangen, Germany 
^Department of Physics, Oklahoma State University, Stillwater, Oklahoma 74078-3072, USA 

(Dated: May 25, 2010) 

We demonstrate that intensity correlations of second order in the fluorescence light of N > 2 single- 
photon emitters may violate locality while the visibility of the signal remains below 1/v^ ?^ 71%. 
For this, we derive a homogeneous Bell- Wigner- type inequality, which can be applied to a broad 
class of experimental setups. We trace the violation of this inequality back to path entanglement 
created by the process of detection. 
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INTRODUCTION 



The demonstration of non- locality for a system of more 
than two particles or even for an EPR state of two parti- 
cles with higher spins has been an outstanding problem 
in the foundations of quantum theory. There are seminal 
papers on the theoretical aspects of this subject, notably 
from Mermin and others [1-5]. These papers pointed out 
the necessity for the use of other Bell- type inequalities [6| 
than the celebrated CHSH or CH74 inequalities [7, ^ in 
case of a system with N > 2 particles or for higher spin 
systems, involving the measurement of N-particle prop- 
erties. However, so far, experimental realizations of these 
inequalities do not seem to exist. 

In this paper we propose as a source for N > 2 par- 
ticles a chain of N independent single photon emitters, 
say trapped ions, trapped neutral atoms, quantum dots, 
or any other equivalent physical system with access to 
similar behavior, in order to test non- locality with more 
than two particles. Using this source, path entangle- 
ment among the emitted photons is created in the pro- 
cess of detection due to the absence of which- way infor- 
mation when registering a photon in the far field of the 
source [y,[lQ|- Employing two point photon-photon corre- 
lations denoted by G^^^(ri, r2) in Glauber's notation [11| 
- a quantity which is easily accessible experimentally and 
has become the workhorse of experimentalists in quan- 
tum information science [12] - we are able to show that 
it is possible to violate a new form of Bell-type inequal- 
ity. In particular, we show that the violation subsists 
for increasing N when the visibility of the photon-photon 
correlation signal continuously reduces to 33%. We note 
that several theoretical and experimental papers have 
dealt with the question of violations of Bell's inequali- 
ties employing the photon-photon correlation function for 
two independent emitters [l3l-[l5j. In this case violation 
of locality has been demonstrated if the visibility of the 
two-photon signal exceeds l/>/2 ^71%, what is recov- 
ered by our results. Our work therefore clearly suggests 
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that without indicating N the magnitude of the visibil- 
ity of the G^^^ -signal can not be taken as a signature of 
non-locality. 



In the following we start to introduce our system of 
N independent single photon emitters. Hereby we as- 
sume that each emitter is initially prepared in an excited 
state which has an appropriate Zeeman degeneracy so 
that each emitter upon spontaneous decay scatters ei- 
ther a right hand or a left hand circularly polarized pho- 
ton. Alternatively, we could consider N two-level atoms 
with a A/4 wave plate positioned in front turning the 
polarization of the scattered photon into the desired cir- 
cular polarization. We further assume two spatially sepa- 
rated detectors in the far field of the source each capable 
of measuring single photon events. The two detectors 
are used to measure G^^^(ri,r2) in order to character- 
ize the two point correlations of the photons emitted by 
our source. As will be shown below, for N > 2 emitters 
we need new Bell's inequalities to reveal the non-local 
behavior of the photon correlations. We derive these in- 
equalities and obtain conditions on the spatial locations 
of the detectors in order to prove the non-local character 
of the photon correlations. These arise from the path en- 
tanglement created a posteriori by the selection of modes 
due to the process of detection [10] . This is a novel aspect 
of our scheme as we do not need to start with a source 
producing entanglement ab initio among the N emitters. 



The paper is organized as follows: in Sec. [Ill we in- 
troduce our light source of N uncorrelated single-photon 
emitters and explain how to describe a joint detection 
measurement of two photons in the far-field of this source. 
In Sec. nil) we recapitulate the well-known set of CH74 
inequalities [3, [8| and explain how these can be violated 
by the probability of finding two photons at two posi- 
tions scattered by our source with N = 2 emitters and 
cannot be violated in case of N > 2. In Sec. IIV( we de- 
rive a new Bell-type inequality [13, Usl which allows to 
prove the non-local character of the correlations among 
the photons scattered by our source for any N > 2. In 
Sec.|V]we finally conclude. 



II. DESCRIPTION OF THE PHYSICAL 

SYSTEM 

A. Setup of N single photon emitters 



a successful detection event of a photon at the detector 
at Tj {j = 1,2), after having passed a polarization filter 
oriented along r]j = sini^jO-+ + cosi^jtr", can be writ- 
ten in case of an even number of emitters in the initial 



We consider the setup shown in Fig. [TJ N single- 
photon emitters regularly arranged in a row at positions 
Ri, R2, . . . , Rat serve as a source for N photons. The in- 
ternal level scheme of the emitters is assumed to be char- 
acterized by a ^-configuration, e.g., Zeeman sub- levels 
with two excited states |e, — 1) and |e,+l), which both 
decay to a common ground state |^,0), accompanied by 
the emission of a a^ or a~ polarized photon, respec- 
tively. For the sake of simplicity we suppose that both 
transitions are equally probable. We further assume that 
for an even number of emitters, the first N/2 atoms are 
initially in the state |e, — 1) and the remaining N/2 atoms 
in the state |e,+l). The initial state of the system can 
thus be written in the form 
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\i^i)=Y[\e.-^)n^ n l^'+l)n. 



(1) 



1: 


|e,-l>^|9,0> 
|e,-l)^|5,0> 


N/2: 


|e,-l>^|g,0) 


N/2+1 


e,+l>^|9,0> 


N-1: 
N: 


e,+l)^|g,0> 
e,+l)^|9,0> 




trap 

FIG. 1: Setup used for measuring the intensity correlation 
function of second order for a source of N single-photon emit- 
ters of known origin and polarization. The photon detectors 
are placed at ri and r2 in the far-field region of the emitters 
and are equipped with polarization filters in front, transmit- 
ting ?7i and 772 polarized light, respectively. 



state ([T]) in the form |1( 



where the subscripts n refers to the atom located at R^. 
For an odd number of emitters, we suppose that the first 
(N — l)/2 emitters are initially in |e, — 1) and the remain- 
ing {N + l)/2 emitters in |e, +1) so that the initial state 
is given by 



N/2 



N 



\^i)= n |e'-l)n^ n l^'+l)n- 



(2) 



Due to the process of spontaneous decay the N three- 
level emitters will scatter exactly ^ {^^^) ^^ ^^^ T 
(^^^^) cr~ polarized photons in the case of an even (odd) 
number of emitters. Alternatively, we could also consider 
N two-level atoms with A/4 wave plates positioned in 
front of each particle which turn the polarization of the 
photons emitted by the atom at R^ forn = 1 . . . ^ ('^^^) 
into t7+ and for n = y ('^^^) • • • A^ into ct~ polarization 
in the case of an even (odd) number of emitters; the only 
prerequisite for our scheme is that a precisely determined 
number of cr^ and C7~ polarized photons of known origin 
is emitted by the setup. 

In order to measure the intensity correlation function 
of second order we locate two detectors at ri and Y2 in the 
the far-field region of the emitters, each equipped with a 
polarization filter in front, oriented along r\\ and 772, re- 
spectively. The operator DN{S{rj),r]j) which describes 



^ n=l 

N 

+ cos^, ^e^^^^|^,0)n(e,+l|) (3) 



where the sum over n takes into account that principally 
each atom could have emitted the recorded photon. Here, 
|^,0)(e,±l| is an atomic operator projecting the atomic 
state |e,±l) onto |^,0), £^0 is the amplitude of the elec- 
tric field, and S{rj) is the optical phase difference between 
photons being emitted by adjacent atoms and registered 
at Tj. In the far-field and using a coordinate system 
where R^ = nRi (with Rq = 0) the optical phase dif- 
ference Sj is given by [16] 



^(r^) :=Sj = kd- 



Ri 



|Ri| 



= kd sin Oj , 



(4) 



where k denotes the wavenumber of the scattered light, 
d the interatomic spacing and Oj the scattering angle as 
shown in Fig.[TJ In analogy, the detection operator for an 
odd number of N emitters acting on the initial state (J2j) 
can be written as 



Dn{Sj,Vj) = -/| 



sin'^. 



n=l 



N 



' \g, 0)„(e, -1| + cost?,- ^ e™^^ \g, 0)„(e, +1| 



(5) 
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FIG. 2: Plot of the intensity correlation function of second order G}^ {6i, —Si, f , f ) (in arbitrary units) as a function of the 
relative phase shift ^i (c.f. Eq. ([7j)). The plot illustrates the signal for different numbers of emitters N (N = 2, 4, 6, 8). 



B. Intensity correlation signal of second order 

With the detection operators Dn{^j^Vj) ^^ hand, we can calculate from Eqs. ([I])-(|5j) the intensity correlation 
function of second order Gj^ ( Ji , (^2 ; ?7i , 772 ) for our system of N single photon emitters. Hereby, we assume in the 
following that the first two out of TV scattered photons are recorded by the two detectors [28]. According to [16] we 
then have 






(6) 



(2) 

As there is a unique correspondence between rfj and dj^ we can write in the following G)y ((^1, (^2; ^1, ''72) also as 
G]y ((^1, (^2; '^1, '^2)- By fixing the orientation of the polarization filters in front of the two detectors identical to 
'^1 = '^2 = f , corresponding to rfj = 1/^/2((t+ + a~) (for j = 1,2), we obtain 



<^N (^1,^2,-, -J 



E^ 



N 



N{N- 



Y.{N-n)cos{n{52-S^)) 



(7) 



n=l 



which holds for even or odd N. This function is illustrated for different N in Fig.[2l In the case that the two polarizers 
are set orthogonal Sit i^i = j and '^2 = ^ (corresponding to rji = l/\^{a~^ + cr~) and 772 = 1/a/2(— cr+ + cr"), 
respectively) we find for an even number of emitters N 
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N/2 



N{N-1) 



- ^{N-2n)cos{n{52-Si)) 



(8) 



n=l 



N/2 N 



N{N -1) 



— ^ ^ {e{N - n - a + 1) cos(n {62 - Si)) e(n - a + 1)) 



Q; = 1 n=l 



where the Heaviside step function Q{x) is defined as 



e{x) :-- 



x< 

1 a;> 



(9) 



In analogy, we find for an odd number of emitters N 
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N-l ^^ 



(10) 



N{N -1) 



— ^ ^ (e(7V - n - a + 1) cos(n {62 - Si)) e(n - a + 1)) 
I 



Q; = 1 n—1 



From Eq. ([7]) we can calculate the visibility Vat 
of the intensity correlation signal of second order 



G^^\5i 



N v^i5^2; f , f ) in case of identically oriented polariz- 



ers. For even or odd N we find 



Vn'--- 



max[G\ ^^ 



^N 



• min[G}, ^^ 



^N 



N 



max[Gjy-'] 



max[G^ ' 



where 

the maximum (minimum) value 



N 



■ rmn[Gjy 
{min[G^^^]) 



37V-4' 



(11) 



corresponds to 
of the function 



G 



(2) 

AT 



G 






^2; f : f )• Eq. (pT]) shows that Vat can be 



uniquely assigned to the number of emitters N . Note 
that Vat represents an ideal theoretical value only; in 
general, experimental uncertainties and insufficiencies 
will influence and decrease the attainable visibility. 

In the derivation of Eqs. ([7|)- (pT]) it has been assumed 
that the two photons measured are the first two pho- 
tons being scattered by our system of N single photon 
emitters. This scenario enables to work with the ini- 
tial states given by Eqs. ([T]) or (J2]). One way to achieve 
this experimentally is to measure all N scattered pho- 
tons and pick out the first two detection events via post- 
selection. Thereby the experimental challenge of measur- 
ing the intensity correlation function of second order for 
A^ possible emitters appears to equal the requirements of 
measurin g th e intensity correlation function of A^th order 
(see, e.g., [16]). We note, however, that our measurement 
scheme requires to resolve the spatial distribution of the 
two-photon correlation signal G]^^ (^1, ^2; '^1, '^2) only and 
that one can make use of a large bucket detector or a lens 
system to detect the remaining TV — 2 photons what sim- 
plifies the requirements. 

Besides experimental challenges, the restriction of de- 
tecting the first two photons bears a major advantage: 
since our system consists of a fixed number of scatter- 
ers, the number of photons contributing to a successful 
measurement cycle is precisely known. Therefore, the in- 
tensity correlation function of NXh order is directly pro- 
portional to the probability of finding N photons. In 
particular, the intensity correlation signal of second order 

(2) 

G]y (^1, J2; '^1, '^2) is related to the detection probability 
of finding jointly the first two photons ^12(^1, ^2; '^1, '^2) 
via 



n2 



-0 ^(2). 



pf2(^i,^2;^i,^2) = ^G)^^((5i,(52;^i,^2), (12) 



where the superscript N denotes the number of emitters 
used in the setup and Cq := ji ^ abbreviates the overall 
success probability to find a single photon at a detector 
with quantum efliciency /i and subtending a solid angle 



III. CH74 INEQUALITIES FOR MULTIPLE 
EMITTERS 



In his seminal paper Bell proved that deterministic lo- 
cal theories with hidden variables are incompatible with 
quantum mechanics [6]. In this Section, we want to apply 
this criterium to investigate whether the photons emit- 
ted by our system of N regularly arranged single photon 
emitters display spatial correlations which are compati- 
ble or incompatible with local deterministic theories. For 
this, we recapitulate a well-known set of homogeneous 
position dependent Bell-type inequalities, the so-called 
CII74 inequalities |3, [sl, [ij , which we then apply for our 
system of N single-photon emitters. 



A. Theory of CH74 inequalities 

Let us denote the continuous set of hidden variables 
by A. The probability of registering one photon out of a 
set of N single photon emitters at a position Yj is then 
determined by j9^(rj. A), where we included the hidden 
variables A in the argument of the single photon detection 
probability p^ {yj). Following the requirement of locality^ 
the joint probability ^^2(^15 ^2, A) of detecting two pho- 
tons at ri and Y2 can be written as the product of the 
two independent single detection probabilities 



^N 



N I 



pf2(ri,r2,A)=p^^(ri,A)-y^(r2,A). 



(13) 



Though A are hidden variables of a deterministic local 
theory and thus unknown, the detection probabilities ob- 
tained when performing a real experiment are determined 
by the ensemble averages over all A 

p^iy,) = j d\g{\)p''{y„\) with i = l,2, 
<2(ri,r2) =ydAfl(A)p^(ri,A)p^(r2,A), (14) 



where g{X) denotes an appropriate weight function of the 
hidden variables. 

Having introduced the single photon and joint detec- 
tion probabilty p^(rj, A) and ^^2(^*5 ^j^ ^)5 respectively, 
the homogeneous CHSH-type inequalities can be derived 
from the following mathematical inequalities [7|, 

- XY <xy-xy' ^x' y^x' y' -Yx' -Xy<{), (15) 

These inequalities hold for any values x^x^y^y^X^Y 
fulfilling < x,x' < X diud < y,y' < Y. Setting 
X = y = 1, so that < x^x'^y^y' < 1, we can then 
identify 

p^(5i,i?i,A) = x, p^{6[,^uX) = x\ p^(^i,oo,A)=X, 
P^(^2,^2,A) = ^, p^((5^,^2,A) = ^^ p^(^2,oo,A) = r.(16) 

where the arguments of the probabilities refer to our 
setup: the jth detector is sensitive to i9j polarized light 



only and is located at Vj (j = 1, 2) where Eq. (j4]) relates 
the detector position Vj to the optical phase Sj. The 
notation 00 indicates that the polarization filter is re- 
moved for the particular measurement. The constraint 
X > x^x' (Y > y, y') is then guaranteed by the so-called 
no- en/iancement condition [3, [sl, [ij* the detection prob- 
ability when using a polarization filter cannot exceed a 
measurement without a polarization filter. Finally, in 
agreement with the requirements of a local hidden vari- 
able (LHV) theory and Eq. (p!3|) . we can write the two- 
photon joint detection probability as 



pf2(<^l,<^2;^l,^2,A)=p^((5i,^l,A)-p^(,52,^2,A). (17) 

Combining Eqs. (p!6|) and (pT|) with Eq. ([15]) we obtain, 
after multiplying the whole expression with ^(A) and in- 
tegrating over A, the following inequality: 



J 



- ^12(^1, ^2; ^1,00) -p^((5i,^2;00,^2)] /Pi2(^l,^2;00,00) < 0. 

I 



(18) 



Hereby, we restricted ourselves to the upper bound of the 
inequalities (fT5]) which allows to normalize the expression 
by an arbitrary function. In the following we choose as 
normalization function the expression ^^2(^17^2; oo, cxo) 
which is a constant, independent of N. This allows in 
particular for a better comparability of the results ob- 
tained in the forthcoming sections. 



Eq. (p!8|) is the position dependent CHSH inequality 
which can be used to investigate the quantum nature of 
the spatial correlations of the photons emitted by our 
source depicted in Fig. [1] (see also [lOj, 114] )• Note that, 
although we are interested in the spatial behavior of the 
two-photon correlation signal, the polarization degrees 
of freedom play a crucial role in the measurements of p^ 
and ^^2- we have to include them necessarily in order 
to satisfy the no-enhancement condition [3, [8|, M]. We 
emphasize, however, that in our investigations we focus 
on the spatial correlations among the emitted photons. 



joint detection probabilities needed in Eq. (fT8|) to 



pUSi,S2;j,1) = f,G%{5^,Sr,H) 







P^SuS^;-,^) = ^G%{5^,5r,-,-) 



q 

^0 



-^G^(<5i,(52;-,— ), 



4'4^ 
4'T 



p\2{5i,52;oo,oo) = Cl 



1 



(19) 
(20) 

(21) 



where we made use of the expressions derived in 
Eqs. (O, (jHj) and ([lO]). Whether or not a violation of 
the position dependent inequality Eq. ([18]) ioi N > 2 
does occur can be verified by inserting Eqs. (fT9|) - (f2T]) 
into Eq. ([T8|) and looking thereafter for the maxima of 
Sn as a function of Ji , 82^ d'^ and 82 . 



B. CH74 inequalities for a system of two 
single-photon emitters (case N=2) 



In order to violate the inequality ([T8]) maximally, it is 
advantageous to adjust the polarization filters such that 
the detection efficiency of the experimental setup is op- 
timized. In the following we thus choose '^1 = '^2 = f 
which yields the best results. With these settings, using 
Eq. ([T2]) and employing the relation ^12(^1, ^2; '^1, co) = 
^12(^1, ^2; '^1, '^2) +^12(^1, ^2; ^1, ^2 + f ) we calculate the 



For the case of A^ = 2 emitters, we find from Eqs. ([7]) 
and ([HJ 



TT TT, E\ 



Gf(5i,52;T,T) = ^(l + cos[J2-5i]) (22) 



4'4 

TT StT, 



8 



G^'n<^i,<52;-,-) = -j-(l-co^[52-5,]) (23) 



so that Eqs. 1^9^ - (|2T]) become 



7iax[S]\[] 



^0 4' 

^0 o- 



Plugging these results into Eq. (p!8|) we obtain 

5'2 = -f C0S[52 - ^l] - C0S[52 ~ ^l] 

+ COS [S2 - S[] + COS [S2 - S[]^ - - 



(24) 
(25) 
(26) 

(27) 



<0. 



Looking for the extrema of So we find the following set 
of parameters (see also |1QLI14|) 



S2 
S2 



5i = - 27r, S'2 - (5i 



•27r, 



(28) 



5[ = -2n, S'2-5[ = l2Tr, 



which lead, in combination with ([27|) . to the following 
inequality with respect to the spatial correlations of the 
photons scattered by two single-photon emitters 



V2- 1 <0. 



(29) 



The inequality Eq. (|29|) is derived assuming an ideal 
visibility of 100% for the two-photon correlation func- 
tions (Eqs. ([22]) and ([23]) ). However, the visibility that 
can be achieved in a real experiment is usually below 
that value due to experimental uncertainties, limited de- 
tector efficiencies etc. Taking a reduced visibility V < 1 

for G^^^((5i,(52; f , f ) and G^^^ ((^1,(^2; f , x) ^^^^ account, 
Eq. dlni) reads: 



5'2 = V2 • V - 1 < 0. 



(30) 



This inequality n iay be violated only if the visibility ex- 
ceeds -^ ^ 71% [io|,[ii,[i3. 



C. CH74 inequalities for a system of multiple 
single-photon emitters (case N > 2) 

For the case of N = 2 emitters, the extrema of S2 
are obtained using the set of analytical expressions for 
^1? ^17 ^2, ^2 provided in Eq. ([28]) . These can also be writ- 
ten in the form 

6i = ai27r, (^2 = (^ + c.2) 2^, (31) 

o 

2 3 

^i = (g + (^3) TT, ^2 = (g + ^4) TT, 

with ai eN {i = 1,...,4). 

In contrast, for A/" > 2, the joint detection probabil- 
ities present in Sn get more involved (c.f. Eqs. ([7|), (|8|) 



-0,1 



^ even N 
^ oddN 



3 



FIG. 3: Maximum values of Sn as a function of the number of 
emitters N: the numerical calculations show that the CHSH- 
type inequality ()18p is violated only for N — 2^ whereas for 
AT > 2 the inequality holds. For AT = 3 and AT = 4 the 
maximum value corresponds to 5'3,4 = 0. 



and ([To]) ). The values for 5i^ 81^82^ 82 giving rise to max- 
ima of Sn were thus determined numerically. This ap- 
proach unveiled that the maxima of Sn for any N > 2 
can be obtained by choosing 



^1 = (^1 27r, 82 = Q^2 27r, 8^ 



Qfa TT, 82 



: Qf4 7r, (32) 



again with a^ G N {i = 1, ..., 4, 0^3, 0^4 ^ 0). 

The results of our numerical calculations for the max- 
ima of Sn (for A^ = 2, ..., 10) are shown in Fig. [3l For 
A/" = 2 we obtain as before S2 = \/2 - I. For A^ = 3, 4 
we find in both cases £'3 = 5*4 = 0. For A^ > 4 the values 
of Sn are displayed in the plot: we see that the behavior 
is slightly different for even N (red stars) and for odd 
N (blue stars). However, we find that a violation of Sn 
appears only for the case N = 2. 

In conclusion, we see from Fig. [3] that Sn cannot be 
violated by the setup shown in Fig. [1] for N > 2 emitters. 
Taking into account the visibility of the intensity corre- 
lation function of second order as derived in Sec. Ill B I (c. f 
Eq. (dH)) this result is in agreement with Eq.(|30|) [10, l^: 
while the joint detection probability for our setup in case 
of N = 2 shows a modulation with a theoretical visibility 
of V2 = 100%, Eq. ([TT]) reveals that Vn drops rapidly 
with N > 2. This is illustrated in Fig. 21 already for 
the case of A^ = 3 the visibility is reduced to V3 = 60%, 
i.e., below the critical value of 1/a/2 ^ 71%. The latter 
was found to be the required value in order to violate the 
Bell- type inequalities ([7, 8], c.f. Eg. ([30]) ). 

Triggered by these results, we will consider a different 
inequality which is more appropriate for our system in 
the next Section. As it turns out this inequality is able 
to prove that the spatial intensity-intensity correlations 
of the photons spontaneously emitted by our source of 
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^ ^ 
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which holds under the condition that < Xi,X2, ^^3 < 1; 
for a proof of this inequahty we refer to [21j (see also 
Appx. [A]). By identifying Xj {j = 1,2,3) again with 
single photon detection probabilities we could speak of 
Eq. (|33]) as an (inhomogeneous) Bell-type inequality since 
it considers both single photon and joint detection prob- 
abilities. However, as motivated in the derivation of the 
CHSH-type inequality above, the experimental require- 
ments can be eased if the inequality under investigation 
involves only detection probabilities of the same order. 
Hence, our goal is to derive a homogeneous Bell-type in- 
equality on the basis of the above Bell-Wigner inequality 
which considers joint detection probabilities only, being 
subject to the same overall success probability. Our pro- 
posal for a new inequality reads 



< Xi X4 — Xi X2 — Xi Xs + X2 X3, 



(34) 



FIG. 4: Plot of the theoretical visibility Vn (c.f. Eq. {[TT])) of 
the intensity correlation signal of second order as a function 
of the number of emitters N for N = 1, . . . , 11; the visibility 
reaches 50% for N = 4 and 33% for Vn -^ 00, respectively. 



and holds for the constraints < xi,X2,X3 < X4 < 1. 
The proof of ([34|) is provided in Appx. [Al 



N single photon emitters are non-local in nature, even in 
the case of N > 2 emitters, i.e., for a visibility Vn < 71%. 



IV. A MORE SUITABLE INEQUALITY FOR 
MULTIPLE EMITTERS 

In the following we introduce a new Bell-type inequal- 
ity which allows to reveal the non-classical nature of the 
spatial intensity-intensity correlations even in the case 
that a two-photon correlation signal with a visibility less 
than 71% is measured. In fact, as will be shown, this 
new Bell-type inequality allows to reveal the non-classical 
character of the two-photon signal even for a visbility ap- 
proaching 33%. To demonstrate this, we start with a dif- 
ferent mathematical inequality based on a so-called Bell 
Wigner-inequality (see, e.g., ^171, |21|). 



In analogy to the foregoing Section, we consider the 
setup with an even (odd) number of emitters N as dis- 
played in Fig. [TJ Again, a photon detection event reg- 
istered at the jth detector is characterized by two pa- 
rameters: the position Vj giving rise to an optical phase 
Sj and the orientation of the jth polarizer rjj which we 
choose to be oriented along l/\^{a~ + a~^) {j = 1,2), 
corresponding to 1^2 = '^1 = f • The latter optimizes the 
overall success of the photon detection probabilities. 

We identify again the parameters of Eq. (|34|) with the 
following detection probabilities 



r.Nf 



p^^((^i,i9i,A) = xi,p^((52,^i,A): 

P^(^3,^2,A)=X3, p^((^4,00,A) = 



= 3^2, 
^4, 



(35) 



A. Derivation of a homogeneous Bell-Wigner 
(HEW) inequality 

The Bell-Wigner inequality can be written in the fol- 
lowing form p/7| 



< xi — xi X2 — xi X3 + X2 xs, 



(33) 



where 00 indicates once more that the polarization filter 
is removed for the particular measurement. The con- 
straint X4 > X3, X2, xi of the inequality ([34]) i s thus guar- 
anteed by the no-enhancement condition [3, [sl, [l^ • ^^e 
detection probability with a polarization filter cannot ex- 
ceed the measurement without a polarization filter. Fol- 
lowing the usual 



assumptions of an LHV theory, we define the joint detection probability exactly as in Eq. (fT7|) . Using this relation 



together with (|S5)) . the inequahty ^^, after multiplying by g{X) and integrating over A, reads 

- p^2{Si,S3;^u^2) +Pi2{S2,S3;^i,^2)] /pU5i,S2;(x>,(x>) >0. 



(36) 



In analogy to the foregoing Section and to provide a better comparability with the results obtained so far we normalized 
Eq. (|36]) again by the factor ^12(^1, ^2; 00, 00) which is independent of N (c.f. Eq. (|2T]) ). In the following, we refer to 
the inequality (|36|) as homogeneous Bell- Wigner (HBW) ine quality. 



B. Violation of the HBW inequality for a system of 
multiple single-photon emitters 

In this subsection, we will test the HBW inequal- 
ity ([36|) for A/" > 2 single-photon emitters by determining 
the minimum values of T^. For this purpose, employing 
Eqs. ([7]), dHI and ([10]), we make use of the joint detection 
probabilities Eqs. ([19]) - (|2T]) and search for the minima 
of Tat. However, as the analyses get involved and analyt- 
ically intricate, we only provide numerical results, which 
we obtained by scanning through the complete parameter 
space of 61^62^63^ S4 for each N separately. 

The results for the minima of Tn {min[TN]) for N = 
2, ..., 10 are displayed in Fig. O It shows that we have 
min[TN] < for A/" = 2, ..., 10. For TV = 3 we obtain the 
lowest value of 7Tim[T3] ^ —0.254 and for A/" = 10 we have 
min[Tio] ^ —0.118. Even though the values of mm [Tat] 
increase monotonously for A^ > 2, numerical calculations 
indicate that they approach zero only for A/" ^ 00. This 
suggests that a violation of the HBW inequalities (|36|) 
can be obtained for any finite A^. 

Note that for A/" = 2 we obtain min[T2] = -0.125 
which sticks out of the overall behavior. We explain this 
outlier by the fact that T2 depends only on three of the 



min[T]\^] 



-0,2 



-0,3 



-^ even N 
^ oddN 



FIG. 5: Minimal values of Tn (c.f. Eq. ([36])) for AT = 2, ..., 10. 
The plot illustrates a steady violation of the HBW inequality 
Eq. §^. 



four parameters (^1,(^2,^3,^4 since, due to destructive in- 
terference, we have Pi2{^i^^4]'f^i^oo) = ^, i.e., a con- 
stant independent of ^i and S4. In contrast, for Tn with 
N > 2 the term p^(5i, ^4; 'i^i, 00) is not a constant and 
thus can be employed to shift min[TN] towards smaller 
values. 



C. Interrelationsship betAveen violation of the 

HBW inequality and visibility of the two-photon 

correlation signal 

Let us again consider the theoretically attainable 
visibility Vat of the two-photon correlation signal 
0^(61,62; f , f ) (c.f Fig. S]). From Eq. ([U]) we know 
that it is given by Vat = 3^_^ which reaches 50% for 
N = 4 and approaches 33% for A/" ^ 00. At the same 
time we can see from Fig. |5] that the HBW inequality 
remains continuously violated when increasing the num- 
ber of emitters N. This shows that the HBW inequality 
Eq. (|36|) can be violated by an intensity correlation sig- 
nal of second order G)^ {Si^Si^ f ' f ) having a visibility 
of below 71%. In fact, for finite iV, our results show that 
a system of A^ regularly arranged single photon emit- 
ters always displays spatial correlations among the scat- 
tered photons which violate the criterion of locality even 
though the visibility of the two-photon correlation signal 
approaches 33%. This clearly demonstrates that without 
indicating N the magnitude of the visibility of the G^^^- 
signal can not be taken as a signature of non-locality. 



V. CONCLUSION 

In conclusion, we investigated the non-local behavior 
of a system of N > 2 particles, i.e., photons emitted by a 
chain of N independent single photon emitters. Path en- 
tanglement among the emitted photons is created in the 
process of detection due to the absence of which- way in- 
formation when registering a photon in the far field of the 
source. Introducing a new homogenous Bell- Wigner in- 
equality and employing simple photon-photon correlation 
functions which are experimentally easily implementable 
in the laboratory we showed that this inequality can be 
violated for any finite number N even though the visibil- 
ity of the two-photon signal approaches 33% in this case. 
The violation of the homogenous Bell- Wigner inequal- 
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ity unambigiously proves the non-local correlations of 
the emitted particles. In contrast, using the well-known 
CH74 inequalities, it turned out that no such violation 
can be obtained for N > 2. For this a visibility greater 
than 71% is required which cannot be achieved for A/" > 2. 



consisting of products of the form Xi Xj {i^j = 1, 2, 3, 4) 
only. Eq. (\A2\i holds if 1 > X4 > Xi, ^2, X3 > is fulfilled. 
For the proof we consider two cases: 

First, we assume X2 > ^i- In this case we can rewrite 
the inequality (|A2[) as 
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Appendix A: Proof of inequality (|34p 



which is valid since both brackets are positive or zero 
due to the fact that X2 > xi and X4 > X2 (note that 

Xi,X2,X3,X4 > 0). 

Second, we assume xi > X2. Here, we make a further 
case differentiation: let us assume xi > xs and rewrite 
the inequality ()A2p as 



In this appendix we prove the inequality (|34|) which is 
an extension of the Bell-Wigner inequality (|33|) (see, e.g., 
[21.]). The Bell-Wigner inequality usually reads 



< Xi — Xi X2 — Xi X3 + X2 X3, 



(Al) 



which is valid under the condition < Xi,X2,:r3 < 1. 

As explained in Sec. |IVl it is advantageous to use the 
inequality (|M|) which reads 



< Xi X4 — Xi X2 — Xi Xs + X2 X3, 



(A2) 



< Xi (X4 - X2) - X3 {Xi — X2). 



(A4) 



This inequality is valid since the first bracket is bigger 
or equals the second due to X4 > Xi and since Xi > X3. 
In contrast, if we assume x^ > xi, we can rewrite the 
inequality (|A2[) as 



X3>a^i 



Xi (X4 - Xs) - X2 {xi - Xs) > Xi (X4 - X3) > 0,(A5) 

where the last inequality holds due to X4 > x^. 
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If other photons than the first two are measured the ex- 

(2) 
pressions for G}^ have to be multiplied by an overall 

factor. This does not change the visibility of the expres- 
sions but introduces different normalization factors when 
considering successive measurements. More information 
regarding this topic is given at the end of the current 
Section. 



